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We study the many-body electronic state created by a time-dependent drive of a mesoscopic 
contact. The many-body state is expressed manifestly in terms of single-electron and electron-hole 
quasiparticle excitations with the amplitudes and probabilities of creation which depend on the 
details of the applied voltage. We experimentally probe the time dependence of the constituent 
electronic states by using an analog of the optical Hong-Ou-Mandel correlation experiment where 
electrons emitted from the terminals with a relative time delay collide at the contact. The electron 
wave packet overlap is directly related to the current noise power in the contact. We have confirmed 
the time dependence of the electronic states predicted theoretically by measurements of the current 
noise power in a tunnel junction under harmonic excitation. 
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Recent years have seen a tremendous experimental 
and theoretical progress in the emerging field of electron 
quantum opticsFollowing the example of optics, the 
quantum nature of electronic transport has been demon¬ 
strated in electronic Mach-Zehnder interferometer— and 
Hanbury Brown-Twiss^”— and Hong-Ou-Mandel^“— inten¬ 
sity correlation experiments. Although quantum optics 
with electrons is in general analogous to the one with 
photons, there are important distinctions between the 
two due to differences in particle statistics, vacuum state 
(Fermi sea vs. photonic vacuum), interaction between 
electrons, decoherence, etc. In particular, a simple con¬ 
stant voltage source can act as a single-electron turnstile^ 
due to the Fermi statistics which is responsible for regu¬ 
lar emission of electrons on a time scale hjeV, where e 
is the electron charge, h is the Planck constant, and V is 
the dc voltage drop over the conductor. 

A step forward towards electron quantum optics has 
been made recently with the realization of on-demand 
electron sources^i^i^iiSrii which can create single- to few- 
particle excitationsThis facilitates the full control 
of the quantum state of electrons in mesoscopic con¬ 
ductors and the dynamical control of elementary exci¬ 
tations using suitably tailored voltage pulsesIn par¬ 
ticular, time-dependent drive creates quasiparticle ex¬ 
citations in the Fermi sea that are single-electron and 
electron-hole pairs whose number and probability of cre¬ 
ation depend on the shape and the amplitude of the ap¬ 
plied voltage^^^ Lorentzian pulses V{t) of a quantized 
area f eV{t)dt/h = N {N is an integer) are special as 
they create N electrons above the Fermi level leaving 
the rest of the Fermi sea unperturbed^i^ Experimentally, 
the presence of electron-hole pairs in the system can be 
seen in the zero-frequency photon-assisted current noise 
power which is increased with respect to the dc noise 
Ievelj2i“i2& More recently, quantum noise oscillations have 
been observed in a driven tunnel junction^i and noise 


spectroscopy using a more complex biharmonic voltage 
drive has been carried out approximating Lorentzian 
pulsesThe creation of single-electron excitations has 
been realized experimentalyi^ and the resulting quantum 
states have been reconstructed using the quantum state 
tomography^i^ 

Even though the progress has been exceptional, a fun¬ 
damental question remains: What is the many-body elec¬ 
tronic state created in the Eermi sea by a voltage drive? 
In this article we find that the many-body state is: 

m = C^'[[(^7T^k + iV^AlB2)\F). ( 1 ) 

k 

Here, = f dSv'^(S)c^S) is the creation operator of a 
single-electron quasiparticle state, Aj^ = f dSu.^ {E)cA£) 
and Bk = J dEu.^{E)c{S) are the operators that create 
electron and hole from the electron-hole pair, (c) are 
the electron creation (annihilation) operators, and |F) is 
the filled Eermi sea. [We have assumed, for simplicity, 
that there is one single-electron quasiparticle created per 
period, which is the case for eV^^/hu: = 1 where Vdc is 
a dc voltage component and u is the frequency of the 
drive.] In addition to a single-electron excitation, there 
is a number of electron-hole pairs created in the system 
(labelled by k = 1,2,...) due to the ac voltage compo¬ 
nent. The probabilities of the electron-hole pair creations 
Pk and the single-electron and electron-hole quasiparticle 
amplitudes y{E) and u^(5) depend on the properties of 
the applied voltage. Eor optimal Lorentzian drive, there 
are no electron-hole pairs created {pk = 0) and the state 
has only single-electron excitations, as expected. 

Quasiparticle amplitudes v and u^ in Eq. © also give 
the time-dependent probabilities of single-electron and 
electron-hole pair creations, |v(t)p and |u^(t)p. The 
time dependence of the wave functions that constitute 
the many-body state in Eq. can be probed by an elec- 
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tronic analog of the optical Hong-Ou-Mandel correlation 
experiment^ where electrons emitted from two termi¬ 
nals with a relative time delay collide at the contact 
When the wave packets arrive at the contact simultane¬ 
ously, their transmission in the same output channel is 
blocked by the Pauli principle, which suppresses the cur¬ 
rent fluctuations. The magnitude of the noise suppres¬ 
sion is proportional to the wave-packet overlap at the 
contact. In the present paper, we have performed this 
experiment and measured current noise power in a tun¬ 
nel junction driven by harmonic time-dependent voltage. 
We have found that the correlation noise as a function of 
a time delay is in agreement with the theoretically pre¬ 
dicted quasiparticle amplitudes of electrons and electron- 
hole pairs in Eq. O- 

Let us consider a generic quantum contact with spin- 
degenerate transmission eigenvalues Tn that are inde¬ 
pendent of energy. The contact is driven by a voltage 
V{t) = Vdc + 14c (^) 5 where Vdc is a constant dc off¬ 
set and 14c(^) is a periodic ac voltage component with 
zero average and the period T = 27r/ci;. The cumulant 
generating function of the charge transfer statistics is 
given hy^ S{x) = 2 E„ TVln[l + /^(l - fR)Tn{e^^ - 
1) + (1 - fL)fRTn{e - 1)]. Here, f^R) are general- 
ized electronic distribution functions in the left (right) 
terminal which depend on two time or energy argu- 
ments; h = f = /r = 

— Here, = fg dt'eVac{t')/h (hereafter h = 1), 
f{£) = 1 - e{£), and fvdd^) = f{£ “ eVdc), where 
0 is the step function. Function / couples only ener¬ 
gies which differ by an integer multiple of cj, f{S\S") = 
T,k afc«fe+m-n/(^'“^‘^“eVdc), where £' = e + noj, £" = 
e + 77110 (0 < e < to), and a„ = (1/T) fg 
In particular, diagonal components of / are given by 
f(.£) = I2^=-oo l«nP/(^-na;-eVdc)- Function/(£:',f") 
is the generalized distribution function of a driven quan¬ 
tum contact in energy representation, while f{£) is the 
stationary nonequilibrium electronic distribution which 
is realized in the junction due to time-dependent drive. 

Next we obtain a decomposition of / into single- 
electron and electron-hole states. As shown in Ref. 
the notion of single-electron and electron-hole pair exci¬ 
tations is related to the eigenproblem of {h, h} = hh+hh, 
where / = (1 — h)/2 and / = (1 — h)/2. For inte¬ 
ger eVdc/^ = there is an Wdimensional subspace of 
{h, h} spanned by N special vectors that are eigenvectors 
of both h and h: hv = — v, hv = v. This subspace corre¬ 
sponds to N electrons injected to the contact per voltage 
cycle. In addition, the operator {h, h} has a series of 
two-dimensional subspaces that are spanned by vectors 
Vq, and v_Q, = /iVq, which are given by hhwo^ = 

The spaces span(vc^, v_q,^) correspond to the electron- 
hole pairs (labelled by k = 1,2,...) created per voltage 
cycle with the probabilities pk = sin^((a/c/2). 

At zero temperature, / and / commute with {h, h}, 
that is, they reduce in the single-electron and electron- 
hole pair subspaces of {h, h}. For simplicity, we restrict 


our consideration to the case A" = 1 in which there is 
only one electron injected per voltage cycle. The eigen¬ 
problem of {h,h} defines a resolution of the identity 
|v)(v| + iA = 1, where |v)(v| is the projector on the 
single-electron state and A = \^ak) (va^ | + |)(v-a^ | 
are the projectors on the electron-hole subspaces. This 
defines a decomposition of / into single-electron and 
electron-hole contributions, 

/= iv)(v|y]/fe, (2) 

k 

where fk = fPk- Similarly, for the completely filled 
Fermi sea of the right lead we have / = fk where 
fk = fPk • The single-electron state v is given by / v = v 
and /v = 0. 

The first term |v)(v| in Eq. (|2]) describes a single¬ 
electron state injected to the contact while fk describe 
the electron-hole pairs. By taking diagonal in time com¬ 
ponents of |v)(v| or fk we gain information on the time 
dependence of single-electron and electron-hole pair wave 
functions. Similarly, by taking diagonal in energy com¬ 
ponents we gain information on the single-electron and 
electron-hole pair contributions in the overall distribution 
function. Indeed, for Lorentzian pulses 14or(0 that carry 
a single charge quantum per cycle, there are no additional 
electron-hole excitations and the time-dependent proba¬ 
bility of the single-electron injection is proportional to 
the drive, |v(t)p = el4or(t)/ci;. This is no longer true for 
a general time-dependent drive where |v(t)p ^ eV{t)luj 
due to the presence of electron-hole pairs. 

Before we proceed with the specific examples, let us 
bring fk in a more transparent form. States Vq, and 
v_Q, in general possess both positive and negative en¬ 
ergy components. We can make a rotation of the basis 
u± = (vq,±v_q,)/v^ in the subspace span(vQ,, v_a) such 
that new basis vectors u+ and u_ possess non-zero com¬ 
ponents only for f > 0 and 5 < 0, respectively. Using 
hva = it is straightforward to check that /u+ = 0 
and /u_ = u_ , as required. The projector Pk in the new 
basis reads Pk = |u^)(u^| + |u^)(u^|. We obtain 

fk =Pfc|u+)(u+| +gfc|u^)(u^i| - iy'pkQk |u+)(u^| + h.c. 

(3) 

and fk = |u^)(u^| {pk = 1 — p/c). Decomposition of 
/ in Eqs. © and m enables us to identify the many- 
body electronic state in the left lead created by the drive. 
We find that f{£',£") = {£')c{£")\^) where |T) 

is given in Eq. O- The right lead is assumed to be 
unperturbed and serves as a reference. 

The physical meaning of the amplitudes u± is mani¬ 
fest in Eq. m and can further be elaborated by taking 
diagonal in energy components of fk. For f > 0 we find 
that fk{£) = P/c|u^(^)P- For f < 0 it is more conve¬ 
nient to consider the distribution of holes f^^ = Pk — fk' 

fk^\^) — Pkp-{^)\‘^• Therefore, u+ and u_ describe 
electrons and holes generated in the system. This no¬ 
tion is also supported in time domain with diagonal in 
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time components fk{t) = |u^(t)p and f^\t) = |u^(t)p. 
Apart from electron and hole states on the diagonal, 
fk contains a lso the off-diagonal terms proportional to 
responsible for mixing of the two, 
see Eq. ©• Electron-hole pairs with pk ^ 0 give no con¬ 
tribution to the transport. On the other hand, for p/c ~ 1 
the electron and the hole from a pair are practically de¬ 
coupled from each other (off-diagonal mixing terms van¬ 
ish), cf. Eq. ([T]). This can also be seen in the cumulant 
generating function which becomes a sum of independent 
electron and hole contributions 

Next we study single-electron and electron-hole pair 
states for different voltage drives, see Eig. [H Let us con¬ 
sider a harmonic drive V{t) = Vdc + Vb cos(co’t) where 
the dc offset A" = 1 is kept fixed while the amplitude 
Vq of the ac component is varied. The drive is char¬ 
acterized by the coefficients an = Jn{eVo/uj) where Jn 
are the Bessel functions of the first kind. In addition 
to a single-electron state, there are also electron-hole 
pairs created in the system and they become more rel¬ 
evant for transport as the amplitude Vq is increased. 
Eor eVo/co < 2 there is only one electron-hole pair with 
Pi = (l/2)(X]^_oo 1^ + -^1 \°'n? - N) in addition to 
the single-electron state injected. Eor eV^juj 1, the 
probability pi practically vanishes and only a single¬ 
electron state remains. Time dependence of the single¬ 
electron wave packet |v(t)p is shown in Eigs. [T]^a-c) in 
comparison to the voltage drive V{t). We find that 
for ac drive amplitudes much smaller than dc voltage 
bias, the electron-hole pair creation is not effective and 
the single-electron wave packet coincides with the drive, 
|v(t)p ^ eV{t)luj, see Eig.^a). Eor ac drive amplitudes 
comparable to or larger than dc offset, the single-electron 
wave packet differs from the drive, see Eigs. [Hb,c). In 
that case the electron-hole pairs become important and 
their wave packets together with the single-electron one 
ensure I(t) = GV{t), where G = (c^/tt) T^. The 

wave functions |u± (t) p of an electron-hole pair are shown 
in Eig. [TJd) for the drive without dc bias {N = 0) which 
does not create single-electron states. The nonequilib¬ 
rium distribution functions f{£) for the voltage drives at 
hand are shown in Eigs. [T]^e-h) (solid curves) together 
with the approximations (dash-dotted curves) computed 
using the most important single-electron or electron-hole 
pair states in Eigs. [TJa-d). Erom Eigs. [Hb,c,f,g) we 
find that electron-hole pairs give a more significant con¬ 
tribution in time domain than in f{£). Indeed, while 
the single-electron wave packets clearly differ from the 
voltage drive, the distribution function is nevertheless 
to a good accuracy given by a single-electron state, 
f{£) « |v(f)p + ^(~^)- This is because for small pk, 
the electron-hole pair functions fk in Eq. ^ have domi¬ 
nant off-diagonal electron-hole mixing components (pro¬ 
portional to ^/pk) which do not contribute to diagonal 
in energy distribution fk{£)^ while they do contribute to 
diagonal in time distribution. 

The time dependence of the electronic wave functions 
can be accessed experimentally by studying current noise 
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FIG. 1. (a)-(d) Wave functions |v(t)|^ and |u±(t)|^ of elec¬ 
trons and electron-hole pairs created by different voltage 
drives V{t). The corresponding nonequilibrium distribution 
functions are shown in the panels (e)-(h). 


power in a setup where two voltage drives with time shift 
r are applied to the terminals, VL(t) = Vdc + Vb cos(ci;t) 
and VR{t) = VL{t — t). This can be viewed as the elec¬ 
tronic analog of the optical Hong-Ou-Mandel (HOM) 
experiments^ in which electron wave packets emitted 
from the terminals with time delay r collide at the 
contactIn the analog HOM experiment, Vbc is the 
static bias voltage between the input and output ports, 
which is in our case only defined with respect to vir¬ 
tual output terminals. Because of the gauge invari¬ 
ance, our two-terminal setup is formally equivalent to 
the case of the voltage SV{t) = VL{t) — VR{t) applied 
only to the left lead with the right lead unperturbed. 
The current noise power as a function of the time delay r 
then reads S' 2 (r) = ^,"^=-00 l’^l'^n[(2eVo/w) sin(wr/2)] 

where So = (b^w/tt) REn with = 1 - 
To express S2{t) in terms of the overlap of the wave 
packets, we proceed as follows. Erom the cumulant gen¬ 
erating function S we find S 2 (t) oc TtU/l — fR)‘^] = 
Tr[(/L - /)2 + (/fi - /)2 - 2(/l - /)(/fl - /)] where 
we have introduced the distribution function / of the 
unperturbed Eermi sea. Here, the first two terms on 
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FIG. 2. Noise S 2 as a function of a time delay r between 
two harmonic signals applied at the leads: (a) iV = 1 and 
eVo/^ = 1-5, 1, 0.5 (solid lines, top to bottom), (b) left lead: 
N = 1, eV^juj — 0.5; right lead: = 0, eV^juj — 2, (c) left 

lead: = 1, eV^juj — 1; right lead: A^ = 0, eV^juj — 2. 

Approximations for S 2 calculated using the overlaps C{t) of 
the wave functions depicted in Fig. [T]are shown for compar¬ 
ison (dash-dotted lines), (d)-(f) S 2 {t) measured in a tunnel 
junction oX uj— 2^ GHz (symbols). Theoretical results are 
shown for temperatures Te = 0 (solid lines) and Te = O.luj 
(dashed lines). 


the right-hand side give the noise when the voltage is 
applied to one lead only while the other lead is unper¬ 
turbed. Both terms give the same contribution to the 
noise Sl = Sr = Sq \n + N\J^{eVo/uj) indepen¬ 

dent of r. The term —Tr[2(/L — /)(/i? —/)] gives the noise 
suppression due to wave packet overlap at the contact. 
The noise reads 

S2{r)/So = {Sl + Sr)/S o - 2C(r), (4) 


noise computed using the most important single-electron 
wave packet overlap in Eq. The corresponding volt¬ 
age drives and the electron wave functions are shown in 
Figs. [U^a-c). 

So far we have analyzed the single-electron wave pack¬ 
ets. To probe the electron-hole states we can use 
the voltages VL{t) = Vdc + VoLCOs{ujt) and VR{t) = 
Vqr cos[c(;(t — r)] which create single-electron wave pack¬ 
ets at the left lead and electron-hole pairs at the right 
lead. The current noise power in this case is given 
by S' 2 (t) = 1^ +-^1 '^n[e^o('r)/w], where 

Vo{t) = - 2VolVorcos{u)t)Y/‘^. In terms of 

the overlap, the noise is given by Eq. where Sr = 
SoE„\n + N\Jl{eVoL/io), Sr = ^ E„ kPn(eWw), 
and C{r) ~ pi?|(v|u^)p is the overlap between the 
single-electron state v at the left lead and the electron 
part of the most dominant electron-hole pair u'^ gen¬ 
erated in the right lead {pr = 0.630). The noise S 2 {t) 
in this case is shown in Eigs. Ei:b,c) ; the correspond¬ 
ing single-electron and electron-hole wave functions are 
shown in Eigs. ma,b,d). 

To verify this picture, we have measured S 2 {t) in a 
tunnel junction under harmonic excitation with uj/27r = 
20 GHz (see Supplemental material). Experimental re¬ 
sults are shown in Eigs.[2jd-f), in agreement with the cur¬ 
rent noise power obtained theoretically. This proves that 
single-electron and electron-hole excitations in Eq. (HD 
can be created by time-dependent voltage and accessed 
experimentally in a noise correlation experiment which 
measures the overlap of the electronic wave functions. 

In conclusion, we have obtained the many-body elec¬ 
tronic state created by a time dependent drive of a quan¬ 
tum contact in terms of single-electron and electron-hole 
quasiparticle excitations. We have confirmed our theoret¬ 
ical predictions by probing the constituent quasiparticle 
states in a HOM-type experiment on a tunnel junction. 
The knowledge of the many-body state opens a way of en¬ 
gineering the required time profile or energy distribution 
of single-electron and electron-hole excitations. Since the 
electronic state in a conductor determines the electro¬ 
magnetic field it generates^ our work can be used to 
produce non-classical states of electromagnetic field, such 
as squeezed or entangled photonic states that have been 
observed recent ly^^i^ 


where C{r) is the overlap that we compute using Eqs. (j2j) 
and ©• Eor p/c ^ 1, the dominant contribution in C{t) 
is the overlap of the single-electron wave functions v(t) 
and v'(t) = v(t — r) injected from the leads. 


|2 


C'(r)«|(v|v')| 


I 


^ V*(f)v(f-T) 


(5) 


The noise S 2 {r) is shown in Eig. [2ja) together with the 
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Supplemental Material for 

^^Electron and electron-hole quasiparticle states in a driven quantum contact” 


I. EXPERIMENTAL SETUP 

In order to generate a single electron state by apply¬ 
ing a periodic voltage bias [V(t) = Vdc + fo cos(ci;t)] on a 
quantum contact, the experiment has to be performed in 
the quantum regime, that is, hw ^ ksT at low tempera¬ 
ture T ^ 100 mK and high frequency u;/2 tt = 20 GHz. To 
provide a good matching to the coaxial cable and avoid 
reflection of the ac excitation, one use a Rq — 50U Al/Al 
oxide/Al tunnel junction as a quantum contact. Indeed, 
the time dependence of single-electron and electron-hole 
wave functions in a driven quantum contact does not de¬ 
pend either on the number nor on the transmission of 
its conduction channels. The experiment is performed in 
a dilution refrigerator where a 0.1 T perpendicular mag¬ 
netic field is applied to turn the A1 normal. A bias tee, 
sketched in Fig. [3] by an inductor and a capacitor, allows 
to separate the dc bias voltage applied on the junction 
from the high frequency part of the setup. The ac exci¬ 
tation is imposed through a directional coupler and the 
current fluctuations emerging from the sample are ampli¬ 
fied by a low noise cryogenic amplifier (noise temperature 
Tat ^ IK). Current fluctuations are band-pass filtered 
between A/ = 1.2 — 3 GHz and the noise power density 
S 2 integrated over the bandwidth of the filter is measured 
with a broadband square law detectors as a function of dc 
bias voltage for different ac excitations. The derivative 
of the noise 882 /d{eVdc) is measured with an additional 
77 Hz, small voltage modulation and a usual lock-in de¬ 
tection. The base temperature Tph of the dilution refrig¬ 
erator significantly increases with increasing ac voltage. 
It implies an increase of the electronic temperature Te 
which is measured by fitting the data of 882 /d{eVdc) (see 
inset in Fig. ED. 


II. CALIBRATION PROCEDURE 


The current noise power associated to the collision of 
two single-electron wave packets generated by two volt¬ 
age drives VL(t) = +Vb cos(ci;t) and VR{t) = VL{t — r) 

is simply given by the harmonic photon-assisted noise 
with Vac = 2eVb sin(c(;r/2): 



tunnel junction 

FIG. 3. Experimental setup for the measurement of the cur¬ 
rent noise power in a tunnel junction under harmonic excita¬ 
tion. 


S2{eVdc,eVac) = X '^ni^) 

where 8 q = Kj/Rq and is the temperature of elec¬ 
trons. Since we measure the amplified current noise 
power density G {82 + 2 kBTN)Kf and we want to ex¬ 
tract its absolute value 82 ^ we need to calibrate the to¬ 
tal gain G of the measurement setup (including the am¬ 
plification and the cable attenuation), the noise added 
by the amplifier Tn and the coupling eVadkjuJ between 
the sample and the excitation line. This calibration 
has been done in two stages: (i) G Af and are de¬ 
duced from the high voltage dependence of the noise 
power which is simply given by the classical shot noise: 


RojeVdc + nhw) 

th ’ 
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S 2 {eVdc ^ eVac^kBTe) = eVdciRo (see red dash line in 
Fig. 13^). (ii) Figure [5] shows the calculated (a) and mea¬ 
sured (b) seconde derivative of the photon-assisted noise 
d‘^{S 2 /So)/d{eVdc)^. One observes maxima correspond¬ 
ing to maxima of J^{eVac/kw) which are approximatively 

aligned with eVac = (0.79±0.05) xfici;+(1.06±0.01) xeV^c 
(black dash line on Fig. [5]). These remarkable points 
are used to calibrate ac coupling and our data are in 
very good agreement with the theoretical predictions for 
eVdclkiJ ^ 0-5. However, we notice a discrepancy for 
eVdclkiJ < 0-5 (red circle on Fig. [Sb) which reveals a 
smaller coupling at low dc bias voltage. It is attributed to 
non-linearities due to Coulomb blockade effects appear¬ 
ing at low temperature and low bias voltage (see inset 
on Fig. IUd). Indeed, in spite of a small effect on the re¬ 
sistance {SR/Ro ^ 1%) the impedance mismatch can be 
emphasized by interference effects. To properly calibrate 
the ac coupling, we have plotted on Fig. [6] the deriva¬ 
tive of the photon-assisted noise as a function of 

normalized ac bias eVad^uj. The fit gives an ac coupling 
6% lower at low bias voltage {eVdc/kLO < 1) than high 
bias voltage {eVdc/kuo >1). 




eVdcIf^ 


FIG. 4. (a) Measured noise power density S 2 and (b) differ¬ 
ential noise power density dS 2 /d{eVdc) for various levels of 
excitation eVolhw. Arrows correspond to the sense of eVo/Hw 
increasing. Inset: Dynamic resistance measured at 80 and 
155 mK. Non-linearities disappear for Te > 400 mK and the 
resistance of the sample is = 50.4 Q. 


eVdc!^ 

FIG. 5. (a) Galculated and (b) measured second derivative of 
the photon-assisted noise (82 jjd(eVdcd a function of 
normalized dc bias eVddkuj and normalized ac bias eVolhw. 
The black dots correspond to the maximum of j\(eVac/kjjj) 
and J^d^aclkjjS)- 



FIG. 6. Galculated (markers) and fitted (solid lines) deriva¬ 
tive of the photon-assisted noise 3 ( 82 !Sq) jd(eVd 8 ) as a func¬ 
tion of normalized ac bias eVaclduj at different dc bias: 
eVddfy^ = 0.061, 1.04, 2.02 and 2.99. 






























